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Deep ocean acoustics, in the absence of shipping and wildhfe, is driven by surface processes. Best 
understood is the signal generated by non-linear surface wave interactions, the Longuet-Higgins 
mechanism, which dominates from 0.1 to 10 Hz, and may be significant for another octave. For this 
source, the spectral matrix of pressure and vector velocity is derived for points near the bottom 
of a deep ocean resting on an elastic half-space. In the absence of a bottom, the ratios of matrix 
elements are universal constants. Bottom effects vitiate the usual "standing wave approximation," 
but a weaker form of the approximation is shown to hold, and this is used for numerical calculations. 
In the weak standing wave approximation, the ratios of matrix elements are independent of the 
surface wave spectrum, but depend on frequency and the propagation environment. Data from the 
IIawaii-2 Observatory are in excellent accord with the theory for frequencies between 0.1 and 1 Hz, 
less so at higher frequencies. Insensitivity of the spectral ratios to wind, and presumably waves, is 
indeed observed in the data. 
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PACS numbers: (43.30.Nb (noise in water), 43.28.Py (Interaction of fluid motion and sound) 



I. INTRODUCTION 

Interest in deep water ambient noise has its origins 
in the 1930s, beginning with attempts to explain micro- 
seisms, persistent and ubiquitous low frequency ground 
oscillations. Their spectrum peaks near 0.17 Hz, roughly 
twice the characteristic frequency of ocean swell, and 
the amplitude correlates with storminess. The most vi- 
able explanation for this phenomena, based on the work 
of Miche,^ was originally proposed by Longuet-HigginsS 
and extended by Hasselmann^ and Brekhovskikhi^ The 
cause is the non-linear interaction of two oppositely trav- 
eling surface waves. Individually, the pressure attenu- 
ates exponentially with depth. However, wave pairs with 
nearly opposed wave vectors collide to excite an acoustic 
wave at approximately twice the frequency which prop- 
agates downward with negligible attenuation. A review 
of the wave- wave mechanism and its seismic effects, to- 
gether with citations to much of the relevant literature, 
can be found in Kibblewhite and Wu<^ 

Here we give a unified exposition of the acoustic field 
radiated by non-linear surface wave interactions. The 
power spectral density matrix (PSDM), obtained from 
autocorrelations and cross-correlations of pressure and 
vector velocity, is found both for an ocean of infinite 
depth, and an ocean layer over an elastic half-space. It 
is shown that the ratios of elements of the PSDM for 
deep observations of the acoustic signal are nearly in- 
variant, depending on frequency and the propagation 
environment but not on the details of the wave spec- 
trum. While the wave spectrum displays dramatic daily 
variability that correlates with wind speed, the propaga- 
tion environment is essentially constant, aside from sound 
speed variations due to seasonal fiuctuations or internal 
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waves. The variability of the wave spectrum manifests it- 
self in fluctuations in the amplitude of individual PSDM 
elements, but not in their ratios. 

The standing wave approximation is based on the fact 
that the sound speed is much greater than the phase 
velocity of surface gravity or capillary waves. Conse- 
quently, only pairs of surface waves with almost exactly 
equal frequency and almost exactly opposing propaga- 
tion directions (hence almost a standing wave) can gener- 
ate a non-evanescent acoustic wave, propagating without 
vertical attenuation (neglecting viscosity). For acoustic 
waves which generate microseisms, the constraint is even 
more stringent; not only must the acoustic wave prop- 
agate without vertical attenuation, it must propagate 
nearly vertically such that the horizontal component of 
the acoustic wave-number is sufficiently small to excite 
microseisms, which have characteristic wave-length ^15 
km. 

It is shown that the standing wave approximation im- 
plies that the ratios of elements of the PSDM are univer- 
sal constants. In fact, this is far too strong a result, and 
we do not expect its consequences to be observed, ex- 
cept perhaps for sufficiently lossy or rough bottoms and 
higher wave frequencies. Because bottom effects negate 
the standing wave approximation, the radiated acoustic 
signal becomes sensitive to very slight deviations from 
the standing wave configuration. 

We show that a weaker form of the standing wave ap- 
proximation will hold since the acoustic source term, un- 
like the acoustic propagation term, is not particularly 
sensitive to variations of the horizontal wave vector fc/i, 
the sum of surface wave vectors, in the region for which 
the acoustic wave propagates without vertical attenua- 
tion, |fc/i| < 27r//c (/ is the acoustic frequency and c is 
the sound speed) . While ratios of PSDM elements are not 
universal constants in the "weak standing wave approxi- 
mation," they will display little variability with changing 
surface conditions, depending only on frequency and the 
propagation environment but not on the details of the 
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surface wave spectrum. 

The theory is verified with data from the hydrophone 
and a pair of three-component seismometers of the 
Hawaii-2 Observatory (H20), located at 5000 m depth 
near 28N, 142W.^ Spectra and cross-spectra for all three- 
hour windows in a 50-day interval of year 2000 have 
been computed. Surface winds, according to the Eu- 
ropean Centre for Medium-Range Weather Forecasts 
(ECMWF), ranged between approximately 2 and 10 
ms^^ during the period, and the significant wave height 
between 1.3 and 2.7 m. There is, indeed, little variation 
in the ratios of the PSDM elements compared with the 
variations in the individual matrix elements; the separa- 
tion between spectra, in dB, is essentially constant. The 
best evidence is from coherency between pressure and 
vertical velocity, which is especially high for frequencies 
less than 1 Hz. The diminution above 1 Hz is attributed 
to bottom effects not embraced by the simple model. 



where Xh = {x, y), q — {qx, qy), q = \q \ and the dispersion 
relation is 



= .99 1 
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(3) 



with surface tension T = .074 N m ^. qgc is the wave 
number of the gravity-capillary transition. 

At next order in the e expansion, the Navier-Stokes 
equations yield an acoustic wave equation and a surface 
boundary condition, with source terms dependent on the 
first-order solution; 



1^ 

C2 dt^ 



9^ + ' 



(4a) 
(4b) 



II. ACOUSTIC RADIATION FROM NON-LINEAR WAVE 
INTERACTIONS 

Starting with Longuet-Higginsj? the most common 
derivation of the acoustic signal radiated by wave-wave 
interactions involves a perturbative solution of the hy- 
drodynamic equations. The expansion parameter char- 
acterizing the non-linearity is proportional to wave slope. 
The leading term is a superposition of plane surface grav- 
ity waves and the second-order term is an acoustic plane 
wave. We give an immensely abbreviated review of this 
approach, closely following Kibblewhite and Wu^ but 
extending the analysis with a unified treatment both of 
pressure and vector velocity and simplifying it by con- 
sidering the bottom as an elastic half-space. Farrell and 
Munk (Ref. 7, Appendix A) give a brief review and ra- 
tionalization of the numerous solutions obtained since 
Longuet-Higgins. 



A. Perturbation equations 

The relevant perturbative solution of the irrotational 
hydrodynamic equations is expressed in terms of a veloc- 
ity potential (j){x, t) and a surface displacement C(a;, y, t). 
The potential and displacement are expanded as 

= 601 -I- 6^02 + • • • (la) 

C = eCl + ■ • • (lb) 

The expansion parameter 6 may be set to 1 at the end of 
the calculation. 

The first-order solution is taken to be an incompress- 
ible fiow corresponding to a superposition of surface plane 
waves; 

0i(f,i) = 

——r / dq Ci{q) exp(i(7 • Xh + qz- ia{q)t) (2a) 

(27r)^ J q 

Ci{x, y,t) = 

1 



dq Ciiq) exp{iq- xh - i(j{q)t), 



(2b) 



where c is the (constant) sound speed. 



B. Acoustic source 

In terms of the Fourier transform of the second-order 
acoustic potential. 



^2i^, kh,z)^ dtdxh(l>2{t,Xh, z)e 



equations (Pa| and (|4b|) become 



+ ^ - kl ] 02 = S{u}, kh, z) 



dz^ 
d uj' 
dz g 
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S{uj,kh,0) 



(5) 

(6a) 
(6b) 



where the source term is [cf. Ref. '6', (4.15)] 

S{uj,kh,z) = 



dtdx, ( ^1 ((V0i)2 



^{-iujt+ikh-Xh) 



Using (|2a]) . 
S{uj,kh,z) 

, . dqdq'a{q)a{q')(l^'^^]C,{q)Uq') 
27rc^ ./ V qq' ' 



(7) 



S{uj - aiq) - a{q'))6\kh - q- q')eM+^"^' 



(8) 



The q' integral is immediately evaluated because of the 
two-dimensional delta function. This gives 



(9) 



S{uj,kh,z) = 
^,jdqaiq).i\k,-q\) (l " ^^f^) 

Ci(g)Ci(4 - - <^{q) - ^(|4 - <f |))e(«+i^^"-'i)^ . 
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C. Solving for the acoustic velocity potential 

In terms of a Green's function satisfying 

j2 , ,2 



where the source S{uj,kh,0) is given in 
The wave elevation spectrum is 



+ ^ - kl) G.,,^(z,.') = -4.<5(z - z') (10a) df.(C(0)r(f))e--, (19) 



d 

dz g 



= 0, 



(10b) 



(cf. Ref. i, Eqns. 4.30, 4.31) the solution of §^ and 
(|6bp is, via Green's theorem, 

1 c2 ~ 

(l)2{uj,kh,z) = S'(a;,fc,j,0)G (^,0) 

47r 5 ' 



_ / dz'^(c.,A:,,z')G^,S,(^,^')- (11) 

The second (bulk) term on the right-hand side of (|lip is 
smaller than the first (surface) term by a factor of order 
gL/<?, where L is a the characteristic depth scale of the 
surface waves, and may therefore be neglected, giving 



^2(w, kh, z) 



±£1 

47r g 



5(w,fc,,0)G^r (z,0). 



(12) 



III. THE POWER SPECTRAL DENSITY MATRIX 

We introduce the 4-vector 

P 

{vQ,vi,V2,v:i) = { — ,v^,Vy,v^), (13) 
pc 

with P pressure and velocity components , Vy in the 
horizontal, in the vertical. The PSDM is defined by 

Af^,„(c^) = ^ [ dT Mt)v,{t + r)) , (14) 

where () denotes averaging over time. The 10 indepen- 
dent elements of ([HI) are evaluated by means of the far- 
field velocity potential 02- With the notation 

{X°,X\X^,X^)^ict,x,y,z), (15) 

one can write 



d 



dXf' 



(16) 



Evaluation of ((T4)) is thus a slight extension of the calcu- 
lation which yields the spectrum of the velocity potential, 

1 f°° 

P^2^7r dTe-'^^{MO,Sh,z)cp;{T,Xh,z)), (17) 

normalized such that di;jF^^{uj) — {(f)\). 
Using one obtains 



b2{(),Xh, z)(i)*2{T,Xh, Z)) 



1 1 

(2^16^ 



dujdkhduj'dk'^ -^{SiLu, kh, 0)S*{io', klO)) 



Giuj,k,,,z)G*{u;',klz)e''^\ 



again normalized such that J dqF^{q) = (C^). Equiva- 
lently, (C(g)C('r')) = {27r)^ F^{q)S^ {q ^ q'). 

Introducing the source spectrum ([9]) into (fT8|) . us- 
ing the elevation spectrum (ITO)) . and assuming Gaussian 
statistics, we get 

{(p2{0, Xh, z)(t)2{T,Xh, Z)) 



dwdfc/j — |G(tJ,fc/i,2;)|^ 



167r2 



\ q\kh-q\ ) 

2Fc{q)F^{kh - q)5{uj - a{q) - a{\kh - q\)y^^ ■ (20) 

Note that one must be careful not to omit a combinatoric 
factor of 2, which arises because 

{amkh~q)Ciq')C{k',-q')) 

= (c(9)r(r))(c(4-<f)C(^;-r)) 

+ {aQ)C{k'H - <7"))(C(4 - q)C{q')) + ■■■ 

The remaining terms, indicated by ■ • • , do not contribute. 
We introduce the function 



S(w, kh) 
87r2 .g2 



dqa{qfa{kh ~ qf 1 



q- {kh - q) 



q\kh -q\ J 

F^{q)F^{kh - q)5[uo- a{q) - aQk^ - q\)) , (21) 



which depends on the surface wave statistics and the dis- 
persion relation. With this definition, the autocorrelation 
of the second-order velocity potential can be written 



2{0,Xh,z)<p^{T,Xh,z)) 

dwd4|G(c^,4,^)|'S(c^,4)e'"" 
duje"^''F^^ . 



(22a) 
(22b) 



Because of the second equivalence, which is the inverse 
of pT]) . it immediately follows that 



F^,{uj)^ / dkh\G{uj,kh,z)f^{uj,kh). (23) 



(18) 



Taking the appropriate derivatives, via equation (|16p . 
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we get 



Mij= j dkhkh,^kh^J\G^J:^{z,Q)\'T.{uJ,kh) 
Mo,r = -- I dkhkhAG^t (z,0)pl](w,4) 



M, 



d 



i,3=i J d-khkh^iG^^^^{z,0) — G*^j^Jz,0)j:{u;,kh) 

Mo.3 = -*^ / dkG^j:Jz,Q)^^Glj:^{z,0)J:{uj,k,) 

The following definitions lead to a more compact no- 
tation 



Ji = |G^,,Jz,0)pE(^,fc, 

^2 = ,-^(z,0)^G;^,^(z,0)E(..,4) (24b) 

— G (z,0)— ^ , 
J4 = 3i(4 = 0) 



(24a) 



^3 = -Gr {z, 0)-G: (z, 0)I](^, fc,0 (24c) 



(24d) 



Moj = -- / dkhkh,j2i [j = f,2) 



-/W^i,3 = « y rffc/i kh,i II2 (j = f , 2) 
A/3,3 = / dkh 1z 



With these, the PSDM may be written 

Ma.o = ^ / dkhli (25a) 
= / d4 (z, J = f , f ; f , 2; 2, 2) (25b) 

CO ... . . _ ^25c) 

(25d) 
(25e) 
(25f) 



IV. THE STANDING WAVE APPROXIMATION 

The collision of two surface waves with horizon- 
tal wave- numbers qi, 52, and corresponding frequen- 
cies cr(gi), cr((72), yields an acoustic wave with fre- 
quency u! — a{qi) + a{q2) and horizontal wave- number 
kh = qi + q2- Thus, the vertical wave- number of the 
acoustic wave is 



M0.3 = / dkh'22 



^=J{u;/c)^-kl 



(26) 



For un-attenuated propagation to the bottom, 7 must be 
real, or fc = < lu/c. This constraint can be re- written 
as 



cr(gi) + g(g2) 
\qi + 92 1 



> c. 



(27) 



Given the disparity between the surface wave phase 
velocity and the much larger acoustic phase velocity, 
(j{q)/q « c, the constraint (P7| implies that qi ~ —q2- 
The allowed deviation from the standing wave-case, 
qi = — 92, expressed as a variation in the relative wave- 
length A A/ A or of the propagation angle A6' of one of the 
opposing surface waves, is of order c^/c, where cq is the 
surface wave phase velocity. Wave spectra are essentially 
constant over such small variations of angle and wave- 
length, which are too small to resolve experimentally. 

Thus, at sufficient depth, the integral over the hori- 
zontal wave number k^ in (|25al - [25^ can be restricted to 
kh < io/c, and the source term Yi{u],kh) can be replaced 
with its value at kh = 0. In the absence of a bottom, the 
squared amplitude of the Green's function |G^ j:^ {z, 0)p 

can also be replaced with its value at kh = for the 
following reason. 

For a bottomless ocean, the Green's function is the 
solution of (jlOap with the surface boundary condition 
(flOb)) . For z < , 



-47rg 



(28) 



(cf. Ref. H, Section 4.2.3). In the region kh < oj/c for 
which there is un-attenuated propagation to the bottom. 



(29) 



All the dependence of |G^_fc^p on kh is in the sub- 
leading term, which is negligible since g/{ujc) is extremely 
small (about 10^"^ at 1 Hz). This is also of order C(^/c, 
since for gravity waves the phase velocity is = g/a 
and, for a standing wave, u) — 2a. Note that the same 
small parameter, g/ujc, determines the dominance of the 
surface term over the bulk term in (jlip . The char- 
acteristic penetration depth of surface gravity waves is 
L = 1/q = g/a'^ — Ag/ut'^, such that the small parameter 
determining the dominance of the surface term over the 
bulk term, gL/c^, is proportional to {g/ujcf' . 

The far-field spectrum is generally computed by re- 
stricting the integration region to kh < to/c and replac- 
ing both source and propagation terms, and G^ , 

respectively, by their values at kh = 0. This is known as 
the standing wave approximation. In light of the above 
discussion, the standing wave approximation is a very 
good approximation in the absence of a bottom, or for a 
perfectly absorbing bottom, with errors proportional to 
the small parameter cq j c. 

To make contact with well-known results, we evalu- 
ate Mo,o in the standing wave approximation. Using J4, 
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Mo,o - ^ / de 



2-K 



dk kli 



TTLO 
C' 
TT W' 

2^ 



— |G„,4(^>0)pS(a.,fc,=0) 

' dqFc{q)Fii-q)Sioj-2aiq)). 



(30) 



The wave spectrum is factored in the usual way, so 
-Fc('?') F(;{q)H{9,q), where / deH{e,q) = 1, and 

/(g) = / d0H{e,q)H{9 + TT,q). (31) 



(32) 



With these definitions, equation ((30)) becomes 
Mo,o(c.) = 7^^F,2(g)/(<z), 



where a; = 2a{q), and w is the group velocity, da/dq. 
The pressure spectrum is then 



(33) 



This result is twice the generally adopted formula (Ref. 
0, Appendix A) , indicating a discrepancy between our ap- 
proach and some previous derivations of the source term. 
The difference, within the spread of prior derivations, has 
no influence on the subsequent analysis of the standing 
wave approximation and the power spectral density ma- 
trix. 

The implications of the standing wave approximation 
upon the elements of the PSDM, other than Mo,o, have 
not been considered previously. 



V. RATIOS OF PSDM ELEMENTS IN THE STANDING 
WAVE APPROXIMATION 



A. Perturbation theory computation 



and all integrals have the limits kh < uj/c. 
The integrals over k^ are trivial, giving 



M 



T^—rJi 



1 §■ 

i 

i 

.1 i 



(35) 



The ratios of PSDM elements are universal constants, 
so we define the ratio matrix. 



R 



M, 



0,0 



2' 

3 


1 

I i 



1 

4 



(36) 



R is independent of both frequency and the ocean wave 
spectrum. We see also that the spectrum of vertical ve- 
locity is 3 dB greater than either of the horizontals. Fur- 
thermore, the normalized pressure is exactly equal to the 
sum of the three velocity spectra on the diagonal. This 
is characteristic of any homogeneous acoustic field in the 
ocean 

Due to rotational invariance, all off-diagonal compo- 
nents vanish except for rg.a and its mirror image. Sen- 
sor calibration determines the accuracy to which the el- 
ements of R can be estimated. A better metric is the 
squared coherency. This will be zero for all off-diagonal 
elements except tq.s. For this term. 



Xo,3 



' 0,3 



The full cross-correlation matrix is 

'^0,0 



'o,i 


' 0,2 


'0,3 




f2,2 


r3,3 




1,2 


1,3 




ri,ir3 3 




r2,2 


'2,3 




''2, 2^3, 3 

?'3,3 - 



(37) 



For a perfectly absorbing bottom and surface wave 
pairs satisfying kh < oj/c, the wave interaction source 
term Ti{uj,kh) and the squared amplitude of the Green's 
function ^^{z,0)\'^ are very well approximated by 

their values at fc/i — 0. Thus, the PSDM elements (|25al - 
l25fD become 



Mo,o = 
Mo,, - 

M^3 
M3,3 

Mo.3 



dkh 

I4 J dkh kh,ikh,j 
--!l4 / dkhkh. 



-1i I dkh kh,i\j ^ - kl 



= 24 



dkh 



u> . 



dkh 



fr2 



(34a) 
(34b) 
(34c) 

(34d) 
(34e) 

(34f) 



B. PSDM elements from incoherent dipoles 

The pressure field from a homogeneous surface layer of 
incoherent and vertically oriented dipoles is 



P{xh,z,t) = duj dxle~'^'F^{x'h)d. 



R 



(38) 



where F is the source amplitude, the subscript h indicates 
horizontal coordinates, and 



(39) 



The more general case, allowing for multiple reflections 
at the surface and bottom, is given by HughesJ^ 
The scaled pressure spectrum, Afo,Oi is 



1 



dS'hV^ix',:, 



R 



(per 



ipc) 



(40) 
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where 'D^{xl^) is obtained from the dipole spectrum, 
{Fui{x)F^{x')) = Duj{x)6uiu>'5{x - x'), and we have 
assumed that we are at sufficient depth z so that 
w/c » 

Since the dipole distribution is assumed to be homoge- 
nous, T)^ is independent of x^. To obtain expressions 
similar to those in (|25al) , we make a change of variables 
from space to wave- number coordinates, -H- fc/i, where 



^£^^{x'^-xkY 



(41) 



is a lucky accident. In general, the Green's function is 
non-analytic at kh = w/c, which is the transition point 
from vertical propagation to attenuation, so there is no 
reason to expect its squared amplitude to be constant in a 
neighborhood of this point. Upon inclusion of a bottom, 
the squared amplitude of the Green's function becomes 
a rapidly varying function on the interval kh = [0,a;/c]. 
The location of the peaks of this function are determined 
by solutions of the eigenvalue equation 



2 



= 0, 



(47) 



Xfi — Xh — z 



Then (l40l) becomes 



A/o,o(w) 



[pc] 



dkh 



(42) 



(43) 



Using P = p(j) and Vi = d4>/dx^, we obtain the rest of the 
PSDM elements: 



Mo,o(w) = 
Mo.iiuj) = 



1 



{per 



c 1 



I dkh 

V,,, I dkh kh,ikh,j 



M3,3(w) 



uj {per 



{pcy 



V,,, i dkhk. 



hKh,: 



V,., I dkh khA\l — ~ k 



{pc) 



I dkh{—- kl 



Mo,3(w) = — 



c 1 



w {pcf 



Vu, / dkh A / ^ - fc: 



(44a) 
(44b) 
(44c) 

(44d) 
(44e) 
(44f) 



and each integral has the limits kh < uj/c 

This agrees precisely with the result p4ap - (j34f|) ob- 
tained using perturbative non-linear wave interaction 
theory in the standing wave approximation, provided 
that the surface dipole spectrum is related to the wave 
spectrum by 



(45) 



VI. BOTTOM EFFECTS 



In general, the effect of the bottom is a modification of 
the Green's function appearing in (|25all25t]) . The stand- 
ing wave approximation which leads to psp is dependent 
on the squared amplitude of the Green's function of the 
operator 



r_ 

dz"^ 



(46) 



being a nearly constant function of kh in the region 
kh < uj/c. The fact that this holds in the bottomless case 



subject to the appropriate top and bottom boundary con- 
ditions. On the real interval kh = [0,aj/c], the Green's 
function has narrow peaks centered about real kh,i cor- 
responding to normal modes and wider peaks centered 
about the real part of complex kh i corresponding to leaky 
modes (see, for example, Ref. [ij. Section 2.1). 



A. Green's function for a layer over a half-space 

We wish to obtain the Green's function G r {z, Zg) 
satisfying 

(^af + ^ - fc^) = S{z ~ zs) , (48) 

subject to the surface boundary condition (for the case 
of gravity waves). 



I G 

9 



(49) 



z=0 



and the mixed Dirichlet-Neuman (impedance) bottom 
boundary condition. 



(G + h{kh)d, 



G 



= 0. 



(50) 



The function h{kh) depends on the characteristics of the 
region below Zbottom = d. 

Assuming constant sound speed c in the water column 
above ^bottom, and taking z^ = (the surface), the solu- 
tion is 



G ,- (z, 0) = 



g (^g-tj(d+z) ^ ^(^^ kh)e'''^'^+^^^ 
(iu2 + ijg)e-'^i'i + R{uj, kh){'jj'^ - 175)6*^"^ 
where 



7 



.2 '^h ' 



and R is the reflection coefficient; 



M + l 



(51) 



(52) 



(53) 



Acoustic radiation and wave- wave interactions 6 



with M = Zb/Z, and the top and bottom layer 
impedances are defined by 



Z = 



Zh 



puj 

7 

ipujh{kh) . 



(54a) 
(54b) 



The function h(kh) entering the bottom boundary 
condition (j50p is obtained from the matching condition 
across the top-bottom interface or across multiple inter- 
faces if the bottom is layered. For a fluid half-space bot- 
tom, one matches the vertical velocity Vz on either side 
of the interface, while for a solid half-space bottom, one 
matches Vz and the components of the stress tensor con- 
taining a vertical index, r^^. This has been done before 
in numerous places (see for example Ref. Il2l ). With- 
out derivation, the result for a fluid layer over an elastic 
half-space is given below. 

For the fluid half-space with density ph and sound 
speed Cb, 



h{kh) = , 
Plb 



lb 




(55) 



For the elastic half-space with density pb, shear wave 
velocity Cs, and compression wave velocity Cp, 



h{kh 



uj^p 



7p 



((2m + A)7- 



4/^fc/i7s 



7p 



— - P 



^2 '^h 



(56) 



where ^ and A are the Lame parameters, related to the 
shear and compression wave velocities by 




(57) 



B. Numerical example 

As a specific example. Green's function for the model 
given in Table |T] was calculated for an acoustic fre- 
quency of 3 Hz (Fig. [T]). A power loss of 6 dB per 
bounce, independent of angle, is assumed for kh > kh,crit 

(fc/i,crit = ^/Cs = 



TABLE L Model used for calculating power spectrum matrix 





Ocean layer 


Half-space 


p(kgm ^) 


1000 


2000 


Cp (ms"^) 


1500 


4400 


Cs (ms~-^) 




2200 


thickness (m) 


5000 


oo 



0.75 



C5 




0.25 



0.000 0.005 0.010 0.015 

Wave number at 3 Hz, radians m'^ 

FIG. 1. The squared amplitude of the Green's function 
Gui,kf^{0, Zhot) on the interval kh ~ [0,cj/c] for a fluid layer 
over an elastic half-space at 3 Hz. Red shows the result for 
a bottomless ocean, which is very nearly flat on this interval. 
The cutoff wavenumber, a;/1500 = .0126, is denoted by +; 
the critical wavenumber, feh^crit = aj/2200 = .0086, with x. 



The figure shows pronounced narrow peaks, 
corresponding to normal modes on the interval 
kh = [kh,ciit,'^/c\ where the reflection coefficient 
has modulus nearly equal to 1 and acoustic energy is 
not lost by propagation through the bottom layer. The 
critical wave number is uj/cg where Cs is the shear wave 
velocity in the half-space. A small power loss per bottom 
bounce is incorporated to give the normal modes kh,i a 
small imaginary component such that these peaks have 
finite height. 

There are also a series of wider peaks on the inter- 
val kh = [0, fc/i^crit], corresponding to leaky modes. A 
general feature of bottom effects is that the acoustic sig- 
nal increases over the bottomless case. Moreover, the 
PSDM ratios can be expected to differ from ([55]) due 
to oscillations of the squared amplitude of the Green's 
function which violate the standing wave approxima- 
tion; it is no longer possible to make the replacement 
Q inside the integrals in ([25a|. 



UJ.kh 



VII. THE WEAK STANDING WAVE APPROXIMATION 

Although the squared amplitude of the Green's func- 
tion is not flat on the interval kh — [0,cj/c] when bot- 
tom effects are considered, the surface wave spectrum 
remains unaltered by bottom effects, except possibly for 
the organ-pipe modes.— With this limitation, one can 
make the replacement 



T.{uj,kh) -> ^iuj,kh = 0), 



(58) 



in p5all25f|) . which we call the "weak standing wave" 
approximation. Since the wave spectrum, E, is evaluated 
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at kh : 
giving 



where 



■^0,0 



"3,3 



■^0,3 



0, it can be taken outside those integrals, 



Mp^i. = S(a;, kh = 0)^^,1. 



T2 J dk,\G^X,iz,0)\' 
dkh kh,ikh,j\G^j:Jz, 0)|^ 



dk 



d 



hkh^G .- (z,0)— G* - (z,0) 



(59) 

(60a) 
(60b) 
(60c) 
(60d) 
(60e) 
(60f) 



and each integral has the limits kh < uj /c. 

The term that depends on the wave spectrum, E, can- 
cels in the ratios of PSDM elements, which therefore de- 
pend solely on S^_y. In the absence of signal other than 
that generated by surface wave interactions, similar argu- 
ments imply that the dependence on the wave spectrum 
also cancels in other ratios of deep ocean acoustic corre- 
lation functions. 

Thus we propose that the PSDM element ratios in deep 
water for the signal generated by wave interactions, while 
not universal constants, are independent of the details of 
the wave spectrum. They will depend solely on frequency 
and the propagation environment, which together deter- 
mine the form of the Green's function G r ■ 



A. Pressure near the bottom 

The effect of the bottom on acoustic pressure is ex- 
pressed as the ratio of the MQfl{z) matrix elements. Us- 
ing the weak standing wave approximation, (|60ap . the 
relative effect of the bottom on pressure is given by 



Mo,o(^) So,o(z) 



(61) 



where ()^ indicates there is a bottom under the model 
ocean. 

Evaluating ([61]) for two depths, we find the solid bot- 
tom raises the pressure at the bottom by an average of 
2.3 dB, and at 500 m off the bottom the increase is 1.6 
dB (Fig. [5]). These values are averages over frequencies 
less than 15 Hz, but the figure shows the smoothed spec- 
trum is flat. Kibblewhite and Wu, using a more detailed 
elastic model, found the bottom lifted the pressure spec- 
trum by 3 dB (Ref. 0, p. 116 and Fig. 7.8). The spectra 
plotted in the figure are equivalent to the function B in- 
troduced by Farrell and Munk to express the effect of the 
bottom on the spectrum of acoustic pressure arising from 
wave- wave interactions on the ocean surface [Ref. i2, (4)]. 




5 10 
Acoustic frequency, Hz 



15 



FIG. 2. The effect of the solid half-space under the ocean 
is a uniform increase in bottom pressure by the factor 1.66 
(black), or 2.3 dB. The pressure increase is slightly less, 1.46, 
500 m above the bottom (red). 



B. Model spectra on the ocean floor 



We have calculated all elements of for the sim- 
ple half-space model (see Table 1). For an observation 
point at the bottom but in the water, the matrix elements 
shown in Fig. [3] are obtained. Smoothing the elements 
over plotted band, 0-15 Hz, gives the ratio matrix. 



R.(zbot) 



-0.1 + 0.4i 
0.23 
. 0.23 
0.29 



(62) 



When smoothed over twice the bandwidth the results are 
essentially the same. If the bottom loss is halved, to 3 
dB, the changes in the matrix are small. Element r3_3 
is reduced by 10%, the other two diagonal elements are 
raised by 10%, and the magnitude of ro.3 becomes 25% 
less. 

The principal effect of the bottom is to reduce the ver- 
tical velocity, relative to pressure, by about 40% with re- 
spect to the case of the bottomless ocean. There are two 
consequences: both ra^s and ro,3 are smaller by about 
this amount. In addition, the coherency between pres- 
sure and vertical velocity is less: Xo,3 = -58, vs. 0.89 for 
the bottomless ocean. The diagonal sum of R is 1.75, 
the difference from 2 reflecting the inhomogeneity of the 
acoustic field near the bottomi^ 



C. Model spectra near the bottom 

As the observation point rises to 500 m above the bot- 
tom, 10% of the water depth, the spectra evolve as shown 
in Fig. m Smoothing the spectra as before gives the fol- 
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FIG. 3. PSDM element ratios r^,,^ — H^,^/Ho,o (black) for a 
point just above the bottom for a 5,000 m water layer over an 
elastic half-space (Table|l]). The short red (blue) line segments 
are the ratios for a bottomless ocean. The red segments termi- 
nate at 3 Hz, the reference frequency of the Green's function 
shown in Fig. [T] 



lowing ratio matrix: 



R(0.9zbot) = 



0.0+.48i 
0.27 
. 0.27 
0.49 



(63) 



The vertical velocity element, r^ ^, has nearly reverted 
to the bottomless value (36), but, surprisingly, the rp.a 
element is 25% less. The diagonal sum is 2.03, reflecting 
the near homogeneity of the acoustic field. 
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FIG. 4. PSDM element ratios (smoothed over 0.3 Hz) at 
the ocean bottom (black, recapitulating Fig. [3| and 500 me- 
ters above the bottom (red). For the diagonal elements and 
frequencies towards the right, the red spectra approach the 
bottomless ocean result, denoted by the red line segments on 
the right margin (panels A, B). The off-diagonal element does 
not (panels C, D). 



VIII. EVIDENCE IN DATA FROM STATION H20 

The theory is applied to bottom data from sta- 
tion H20,— located at 5,000 meters in a thin-sediment 
area of the Pacific, mid-way between Hawaii and Cal- 
ifornia. The instrumentation consists of a buried 
(0.5 m) Guralp seismometer (Guralp System Limited's 
CMG-3) and Geospace geophone (GTC, Inc.'s Geospace 

TM 

Technologies HS-1), and a hydrophone a little off the 
bottom. The two velocity sensors differ in the way the 
motion of the inertial mass is sensed: the Guralp is bet- 
ter at very low frequencies, the Geospace at very high, 
but they are comparable over our analysis band. In ad- 
dition, assimilated surface winds for the H20 location 
were provided by the European Centre for Medium-range 
Weather Forecasts from the ERA-interim (ECMWF Re- 
analysis) results. 



H20 data (hydrophone, channel HDH; Geospace, 
channels EHZ, EHI, EH2; Guralp, channels HHZ, HLl, 
HL2), obtained from the Incorporated Research Institu- 
tions for Seismology (IRIS) Data Management Center, 
have been studied for years 2000-2002, inclusive. Spectra 
were calculated for three-hour windows to a resolution of 
0.1 Hz, giving about 2000 equivalent degrees of freedom. 
The spectra have been examined for the whole interval, 
but this discussion is restricted to days 200-250 in year 
2000. Spectra were calibrated according to the nominal 
transfer functions on file at the IRIS Data Management 
Center. Various small inaccuracies were discovered and 
corrected. For example, there were discrepancies at the 
times the gain was changed. The spectra for the vertical 
component of the Geospace geophone were raised by 3 
dB. This makes them consistent with spectra from the 
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vertical component of the Guralp seismometer, is com- 
patible with the instrument noise model, and corrects a 
demonstrable error in its nominal generator constant. 



A. Bottom acoustics and surface winds 

It is generally observed that bottom acoustics, from 
frequencies less than 1 Hz to frequencies above 30, are 
highly correlated with surface wind.^"^ Typical are the 
profiles of the acoustic power at 0.5 Hz (Fig. 5) and 1.76 
Hz (Fig. 6). The lower frequency was chosen because 
the coherency between pressure and vertical velocity is 
greatest. The higher was selected to be on a relative 
extremum of the coherency and away from any sediment 
resonance. 

The vertical scales in Fig. [SJ Panel A, have been ad- 
justed so that the variation in wind (black) is about as 
large as the variation in acoustics (colored). In Panel A 
of Fig ini the wind has about half the range of the acous- 
tics. Comparing axes, the slope in log acoustic power is 
roughly 2.5dB/(ms~^) at these frequencies. It has pre- 
viously been shown (Ref. 0, Fig 2) that between 6 and 
30 Hz the slope of the spectrum is in the range 2.7 - 2.9 
dB/(ms-i). 

In both figures, the flatness of the components of the 
velocity ratios (Panels B, C) shows that scaling vector 
velocity by normalized pressure is effective at reducing or 
even eliminating the correlation with surface wind. This 
is as expected from the foregoing theory, to the extent 
the wind is a proxy for the waves. 




215 220 225 230 235 240 
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B. Off-diagonal matrix elements 

The off-diagonal elements of the PSDM are the most 
diagnostic of the wave-wave acoustic field because, ex- 
pressed as squared coherency (see Eq. I57)) . they are in- 
dependent of sensor calibration. The element Xo 3 is the 
only off-diagonal element that does not vanish for the 
elementary models we are testing. 

The spectrum of Xo,3 (Fig- 13 shows a profound dis- 
continuity at 1.03 Hz, which is just below the gravest 
sediment resonance (see Fig. 10). The two effects are 
presumably related. For lower frequencies the measured 
coherency is even higher than the model (0.8 vs 0.58, 
see Eqs. 37 and 62). The phase of the cross-spectrum 
(Fig [71 Panel B) in the high-coherency band is about 
—25°, far from the model expectation of 104° (62). How- 
ever, the phase of the cross-spectrum depends on the 
phase response of both instruments, and, as explained in 
the Discussion, there are questions about the hydrophone 
transfer function. 

At higher frequencies the coherency wobbles around 
0.1-0.2, and there are large swings in phase. There is 
pronounced rippling in the coherency, with a periodicity 
of approximately 6 cycles/Hz. These are likely organ 
pipe modes, a possibility previously considered but not 
embraced.— 

In theory, the other five off-diagonal elements vanish; 
in practice, they nearly do so for / < 1 Hz and are not 



Year 2000 day 



FIG. 5. Profile of acoustic power at 0.5 Hz for Guralp data 
with dash red pressure (scaled by (pc)^), and red, blue, green 
denoting the vertical and two horizontal components of veloc- 
ity (A). ECMWF wind (black) is referenced to the right axis. 
Panels B and C are the the three ratios of velocity to scaled 
pressure for the two instruments. Blanks indicate sections 
where there were no data or where the value was dropped on 
a first difference screen. The plus symbols at the bottom of 
each panel denote windows at days 218.125 (U = 3.3) and 
232.125 (U — 10.1) for which detailed results are presented 
below. 



large at higher frequencies (Fig. |8|). The spike in Xi,3 
near 2.5 Hz corresponds to one of the P-SV modes in the 
sediments (see Fig. [TO)) . 



Dividing the velocity spectra by the (scaled) pressure 
(Figs 5 and 6, panels B and C) is effective at obliterating 
the wind correlation. Thus, it is not surprising that the 
cross-spectrum of pressure and vertical velocity is equally 
quiescent (Fig. 9). Another perspective on the residual 
wind signature is visible in comparing the pairs of spectra 
plotted in Figs. 7 and 8. Red is the coherency under a 
strong wind, and black dash under a weak wind. The 
pairs of curves are virtually indistinguishable. 
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FIG. 6. Profiles of acoustic power at 1.76 Hz for Guralp data 
(colored) and ECMWF wind (black) (A). Panels B and C 
are the the three ratios of velocity to scaled pressure for the 
two instruments as in Fig. 5. See its caption for additional 
details. Note that the acoustics in Panel A more closely follow 
the dips in the wind than is the case at the lower frequency 
of Fig. El 



C. Diagonal matrix elements 



The power spectra of the velocity components down 
the diagonal of the PSDM (Fig. ^ are relatively fiat for 
frequencies above 1 Hz but fall precipitously below that. 
Excluding sediment resonances, the spectra are in ac- 
cord with the spot measurements obtained from smooth- 
ing the results at 0.5 Hz and 1.76 Hz displayed in Figs. [5] 
andini respectively, and listed in Table [Hi The steep drop 
below 1 Hz we tentatively attribute to a hydrophone cali- 
bration error (see Discussion). The bumps on the spectra 
of both components at 1 Hz and 2.5 Hz are attributed 
to P-SV sediment resonances, as is the broader peak at 
4.25 Hz for ri,i (blue).ii 




180 




2 3 
Frequency, Hz 



FIG. 7. Pressure and vertical velocity are highly correlated for 
frequencies less than 1 Hz (A), with phase angle of —25° at 0.5 
Hz (B, Table HI)) . The coherency at higher frequencies, though 
less, is still significant, but because of the lower coherency, 
the phase fluctuates more. These Guralp exemplars include 
a high-wind window (red, day 232.125, U=10.1 ms~^) and a 
low-wind window (black dash, day 218.125, U=3.3 ms~^). 



D. Summary 

The model yields a PSDM ratio matrix for which all 
elements are independent of both frequency and wave 
spectrum. This is true whether they are expressed in 
natural units or as squared coherency. 

The measured PSDM ratio matrix is, indeed, insensi- 
tive to wind, which we take as a proxy for waves, but 
strongly depends on frequency. There is a sharp tran- 
sition at 1 Hz. For lower frequencies, the coherency of 
the off-diagonal elements is closer to the model than for 
higher. A sharp transition in the diagonal matrix ele- 
ments occurs at this same critical frequency. Below it, 
the ratio elements plunge by 30 dB, which is entirely 
caused but the sharp rise in the pressure spectrum. Al- 
though we have suggested this is a calibration matter, the 
correspondence in frequency of the two effects suggests 
they are linked in the physics. Above the transition fre- 
quency, the diagonal elements are flatter, which is more 
in accord with the model. The variations are within 5 dB. 
Peaks of up to 10 dB, which are attributed to P-SV res- 
onances in the sediments, are discounted. However, the 
levels are much larger than the model results, an effect 
also attributed to energy trapped in the sediments. 

The important numeral results are highlighted in Ta- 
ble [Hi For six representative matrix elements (column 
1), this shows the model values (columns 2) and the ob- 
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FIG. 8. For the same windows displayed in Fig. [71 the co- 
herency between pressure and horizontal velocity (A) and be- 
tween vertical and horizontal velocity (B) is low everywhere, 
but especially for frequencies less than 1 Hz, where Xo,3 is 
highest (see Table |TI1). 
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FIG. 9. The coherency between pressure and vertical velocity 
(Geospace sensor, three representative frequencies) is virtu- 
ally constant, irrespective of wind speed. For 0.5 Hz, there is 
a slight dimple near day 218 when the wind was low. Results 
for 1.67 Hz around day 235 are unusual, as is the scatter be- 
tween days 240 and 245 at 1.76 Hz. (Gaps are due to missing 
data or spectra dropped on a first difference criterion.) 



served values (columns 3-6) for both sensors and the two 
representative frequencies. 

For the diagonal elements (top rows), the values for 
0.5 Hz are orders of magnitude smaller than the model, 
and the values at 1.76 Hz (and above) are 10 to 40 times 
larger. The off-diagonal elements (bottom three rows) 
fit the model better at 0.5 Hz than at 1.76, although 
in neither case has the result been smoothed over the 
prominent ripples. 

The frequency dependence observed in the coherency 
of the off-diagonal matrix elements is inconsistent with 
the universal constants predicted by the standing wave 
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FIG. 10. Guralp matrix elements ri,i (blue) and rs^s (red) 
for day 232.125 {U = 10.1 ms"^). The adjacent black dash 
curves are the spectra for day 218.125 {U — 3.3ms~^). The 
black bars are averages of the power at 1.76 Hz, exemplified 
in Panel B of Fig. (6] The short bars on the right are matrix 
values from (|62p . Spectra for the Geospace instrument are 
similar. 



TABLE II. Selected elements of the spectral ratio matrix at 
0.5 and 1.76 Hz for both H20 seismic sensors. The 12 values 
for the diagonal elements rij were calculated by averaging 
the spectrum over the 50 days, for which shorter segments 
are plotted in Figs. 5 and 6. The values of the off-diagonal 
elements were calculated by similar smoothing in time at the 
appropriate frequencies (for the Guralp, see Figs. 7 and 8). 





Theory 


0.5 Hz 


1.76 Hz 




(EU 


Guralp 


Geospace 


Guralp 


Geospace 


ri,i 


0.23 


1.86 X 10"^ 


1.39 X 10"^ 


8.29 


6.28 


r2,2 


0.23 


1.54 X 10"^ 


1.57 X 10"^ 


5.82 


6.12 


r3,3 


0.29 


1.86 X 10"^ 


1.70 X 10"^ 


2.21 


2.07 


ro.i 





< .01 


< .01 


< .01 


< .01 


ro,3 


-0.1 + 0.4i 


.81 - .43i 


.82 - .30i 


.32 - A8i 


.32 - .85i 







< .01 


< .01 


< .02 


< .02 



approximation in the case of an ocean layer resting on an 
elastic half-space. Thus, other effects, possibly bottom 
scattering, are influencing the wave-generated sound at 
this site. On the other hand, the insensitivity to wind, 
as predicted by the weak standing wave approximation, 
is strongly upheld by these data. 



IX. DISCUSSION 

We have shown that the theoretical predictions of the 
standing wave approximation are extremely strong, yield- 
ing PSDM ratios which are universal constants. While 
it is possible that these PSDM ratios are observed under 
certain conditions, such as a high loss bottom, they are 
certainly not observed in the H20 data. 
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Bottom effects vitiate the standing wave approxima- 
tion, leaving in place a weaker version of the approxi- 
mation, the predictions of which are observed at H20. In 
particular, the PSDM ratios depend on frequency, but are 
insensitive to changes in the surface conditions. Direct 
comparison of the PSDM ratios to an elastic half-space 
model show less agreement. This is due, presumably, 
to the oversimplification of the bottom and propagation 
conditions in the model, as well as to instrument calibra- 
tion error. 



A. Bottom interaction 

The primitive model of the ocean's bottom replicates 
the essential features of more elaborate approaches. Pres- 
sure is increased a few dB because of the bottom, and 
this increase is uniform over the band 1-30 Hz, at least. 
The augmentation diminishes slightly as the observation 
depth lessens. Given a quantitative estimate for the 
influence of the bottom, overhead wave properties can 
be inferred from (corrected) deep pressure observations 
through application of (33), although there is a factor 2 
discrepancy between this result and some other theories. 
The model of the ocean bottom will need to be extended 
to account for the layer of ocean sediments before velocity 
observations can be used for the same purpose. 

Pressure and vertical velocity are most coherent for 
frequencies below 1 Hz; indeed, the observed coherence 
exceeds the model result. The coherence falls above 1 Hz. 
The sediments overlying the basement were not modeled, 
but the decrease in coherence is attributed to P-SV en- 
ergy trapped in this layer. This same effect can explain 
the weak but significant rise at 1 Hz in the coherence be- 
tween horizontal velocity and both pressure and vertical 
velocity. 

We view the peaks in the velocity spectra as evidence of 
vertically polarized shear waves, naturally excited when 
a pressure wave in an acoustic medium is incident on an 
elastic medium. The extremely low shear velocity of the 
surficial sediments has two effects: the transmitted SV 
rays are nearly perpendicular to the boundary and they 
are strongly polarized in the horizontal direction. The 
resonance occurs when the rays are efficiently reflected 
by internally layering in the sediments. Zeldenrust and 
Stephen^^ applied the theory of Godin and Chapman^^ 
to interpret the resonances at H20 as evidence of a chert 
layer approximately 13 m below the bottom, about in 
the middle of the 30 m of sediments. Scattering is a 
further complication. Any zones in the elastic medium 
with strong impedance contrasts and rough boundaries 
will scatter both SV and SH energy. 



B. Influence of the ocean's sound speed profile 

Taking the ocean to have constant sound speed ignores 
the refraction due to the sound channel (e.g. Ref. 0, Fig. 
3). However, for a surface layer of incoherent dipoles, 
the bottom signal is dominated by the source region with 
diameter six times the water depth fRef [Tsl Eq. 7). 



C. Inference of the wind-wave spectrum 

Bottom acoustic observations are beginning to be used 
to estimate the spectrum of ocean surface waves. To 
do this, within the framework of acoustic radiation from 
wave- wave interactions, two corrections are necessary. 
Allowance must be made for the bottom interaction and 
the overlap integral. 

The correction of the pressure spectrum for bottom 
interaction appears to be straightforward. The bottom 
elevates the pressure a few dB at the bottom, with a 
slight decrease moving up the water column away from 
the bottom. The bottom effect, in fact, is smaller than 
the factor-of-two discrepancy between our derivation and 
some previous results. 

It is just as significant that the bottom correction does 
not depend on frequency. Thus, the theory for the bot- 
tomless ocean can be used to infer the slope of the wave 
spectrum from the slope of the pressure spectrun*^ with 
no further corrections. 

Spectra of bottom velocity at this site are contami- 
nated by sediment effects for / ~ 1 Hz. Horizontal is 
more affected than vertical, and models incorporating 
the sediments will be required to back out this contam- 
ination. Data at lower frequencies may still be applica- 
ble. However, all bottom sensors are equally useful for 
relating changes in the spectrum to changes in overhead 
waves. 



D. Sensor calibration 

1. Hydrophone 

When the nominal transfer function of the H20 HDH 
hydrophone is adopted, spectral ratios of velocity to 
(scaled) pressure, drop precipitously for frequencies less 
than 1 Hz (Fig. [TU|) . At 0.5 Hz, the observed spectral 
ratios for both seismic sensors are more than 100 times 
smaller than theory (Table HI]) . 

There is more evidence that the nominal gain of the hy- 
drophone is too high at low frequencies. Velocity spectra 
at low frequencies have been calculated from ECMWF di- 
rectional wave models, and they are in reasonable agree- 
ment with the H20 seismic observations (Ref. '7, Fig. 
4). Pressure spectra calculated from the same models 
are orders of magnitude less than observations. 

In addition, the (scaled) pressure spectrum at 0.5 Hz 
varies between -95 dB and -105 dB at H20, depending 
on overhead wind (e.g. Figs. [3 |6]). Observations at 
the Aloha Cabled Observatory, scaled similarly, range 
between -113 and -128 dB, some 20 dB lessil 



2. Seismometers 

The similarity between the spectra of all six seismic 
channels at low frequency is in accord with the model 
(cf. panels AA in Figs. [5l IH and Fig. [TO\j and indicates 
relative calibrations accurate to a dB (once spectra of the 
Geospace vertical have been lifted 3 dB). However, the 
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phase of the cross-spectrum between pressure and verti- 
cal velocity differs for the two instruments for frequencies 
less than 1 Hz. Sediment effects preclude applying these 
checks at higher frequencies (e.g. panels BB in Figs. [SI 

El). 
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